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!•  IMEODUCTION^THB  whistler  nozzle 

The  phenomenon  of  entrainment  is  a process  through  which  a 
moving  fluid  increases  its  mass  flow,  by  "grabbing  up"  the 
medium  it  is  flowing  in.  The  moving  fluid  is  called  the 
working  fluid,  and  might  be  the  primary  water  jet  in  an 
ejector  pump.  The  medium  might  be  any  other  fluid,  say  oil, 
in  the  case  of  the  ejector  pump. 

The  entrainment  process  has  found  practical  use  for 
quite  a long  time.  Analytical  descriptions  of  the 
entrainment  process  exist;  however,  none  is  general  enough  to 
describe  all  cases  of  practical  interest.  Even  the 
relatively  simple  case  of  a steady  jet  in  still  air  is 
rather  difficult  because  of  the  many  configurations  which 
one  would  lake  to  describe.  Thus  interest  and  research  in 
this  area  continues,  perhaps  with  more  intensity  than 
previously  due  to  the  recent  interest  in  new  V/STOL 
configurations  which  would  employ  this  process.  One  such  is 
the  Navy’s  Xi’V  - 12A,  which  is  expected  to  derive  a large 
amount  of  thrust  through  the  augmentation,  or  enrrainment, 
of  its  primary  thrust  using  wing  mounted  ejectors. 

The  success  of  concepts  such  as  the  augmenter  wing 
depends  cn  many  factors,  an  important  one  being  the  rate  of 
entrainment  of  the  primary  jet,  and  another  being  its 
efficiency.  This'  rate  of  entrainment  is  the  rate  of 
increase  of  jet  mass  flow  with  axial  distance  away  from  the 
jet  nozzle.  The  rate  of  entrainment  and  nozzle  efficiency 
are  characteristics  which  vary  with  nozzle  configuration. 

As  mentioned  earlier  , there  is  a large  variety  of 
steady  nozzle  configurations.  Lobe,  slot,  and  hyperroixihg 
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nozzles  are  some  which  have  drawn  recent  interest  for  their 
entrainment  capabilities,  one  should  not  neglect  unsteady 
nozzle  configurations,  hov/ever,  as  such  effects  as  swirling 
the  primary  flow,  and  oscillating  the  primary  flow  in  a bi-r 
stable  manner  (fluidic  oscillator),  also  show  promise  in 
their  ability  to  entrain,. 

One  of  the  newest  nozzle  configurations  to  surface  is  an 
unsteady  nozzle.  Perhaps  the  simplest  in  construction  of 
all  unsteady  nozzles,  it  is  illustrated  in  figure  1.  The 
whistler  nozzle,  as  it  is  called,  produces  a loud  pure  tone 
while  in  operation.  Additionally,  it  creates  a relatively 
large  rate  of  entrainmenti  It  is  the  whistler  nozzle  which 
is  the  subject  of  this  thesis. 


I 


II.  BELATED  £AST  WORK 

The  earliest  work  done  on  the  whistler  nozzle  was  by  Hill 
and  Greene  (reference  1) . In  1974,  they  cited  the  discovery 
of  a new  phencmenon,  characteristic  of  a new  device  which 
they  named  the  whistler  nozzle.  Their  basic  configuration, 
as  illustrated  in  figure  1,  consisted  of  a convergent 
section,  followed  by  a section  of  constant  area  followed  by 
a step  change  in  cross  sectional  area  to  another  section  of 
constant  area.  The  length  of  the  Tast  section  could  be 
varied  to  produce  loud  pure  tones  of  varying  frequency, 
through  the  excitation  of  a standing  acoustic  wave  in  the 
final  section,  much  as  in  an  organ  pipe.  Hill  and  Greene 
observed  that  the  whistler  produced  a greatly  increased 
mixing  rate,  and  decided  that  this  was  due  to  acoustic 
stimulation  of  the  jet.  Increased  mixing  rate  of  an 
axisymmetric  jet  due  to  acoustic  stimulation  is  a process 
that  had  been  observed  earlier  by  Crow  and 
Champagne  (ref  erence  5). 

In  addition  to  the  basic  axisymmetric  configuration, 
many  other  configurations  were  found  to  produce  the 
whistling  phenomenon.  A few  are  illustrated  in  figure  2. 
Common  to  all  of  these  was  that  the  step  in  cross  section 
had  to  extend  completely  around  the  jet.  Also  noted  was  the 
apparent  presence  of  a separation  reattachment  cycle, 
occurring  at  the  nozzle  lip; 


Nc  further  study  seems  to  have  been  performed  concerning 
the  whistler  nozzle,  save  a report  in  1975  by  Hill  and 
Jenkins  (reference  2),  in  which  they  reestablished  the 
operating  characteristics  of  the  nozzle,  and  its  ability  to 
entrain  more  air. 
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Hill  and  Greene's  work  has  been  summarized,  since  their 
original  report,  in  comprehensive  reviews  of  the  elector 
state  of  the  art  by  Schum  (reference  3),  and  by  Viets 
(reference  4) . These  reports  call  attention  to  the  whistler 
nozzle's  promise  in  rate  of  entrainment,  and  the  lack  of 
data  regarding  the  whistler  thrust  efficiency. 

Although  no  further  research  directly  involving  the 
whistler  is  available,  much  work  has  been  done  in  the  area 
of  free  jets,  and  entrainment,  which  might  prove  valuable  in 
a study  of  the  whistler  nozzle. 

An  obviously  important  area  of  related  work  is  the 
sensitivity  of  free  jets  to  acoustic  stimulation.  Probaoly 
the  best  known  and  most  definitive  work  in  this  area  v;as 
performed  by  Crow  and  Champagne  (reference  5) . By  mounting 
a speaker  inside  the  plenum  chamber  they  managed  to  excite  a 
free  axisymmetric  jet,  at  various  freguencies.  The  effect 
of  the  sound  was  to  produce  a velocity  fluctuation  in  the 
nozzle  exit  plane.  It  was  observed  that  the  turbulent 
structure  of  the  jet  was  sensitive  to  the  acoustic 
excitation.  Through  the  use  of  high-speed  schlieren 
photography  it  was  observed  that  under  acoustic  stimulation 
the  turbulent  structure  of  the  jet  developed  a large  - scale 
vortex  structure,  as  contrasted  with  the  small  scale  vortex 
sheet  surrounding  a normal  jet.  Attendant  to  the  formation 
of  the  large-scale  vortices,  an  increase  in  the  mixing  rate 
was  also  observed. 

Others,  such  as  Bevilagua  and  Lykoudis  (reference  6) , 
have  observed  the  turbulent  structure  in  steady,  free  jets. 
The  presence  of  a large-scale  vortex  structure  has  been 
observed  in  such  jets  at  low  Reynolds  numbers. - It  appears 
that  the  large-scale  structure  dissappears  as  Reynolds 
number  increases.  However,  according  to  the  findings  of 


Crow  and  -Charapa9ne,-  it  may  be  reexcited  through  acoustic 
stimulation.  It.  may.  be  con'jectured  from  experimental 
evidence  that  the  entrainment  process  is  a combination  of  a 
large-scale  “scooping  up!‘  of  the  ambient  fluid,  and  a small- 
scale  nibbling  away  at  the  ambient  fluid.  These  entrainment 
processes  wculd  correspond  to  the  large  and  small  scale 
vortex  structures  ar.ound  a free  jet. 

A further  area  of  related  research,  perhaps  not  obvious, 
is  the  well  known  Coanda  effect.  The  ’ basic  whistler 
configuration  sets  up  a flow  situation  somewhat .akin  to  that 
of  a bistable  fluidic  element,  or  fluidic  oscillator.  The 
difference  between  the  two  cases  is  that  the  fluidic 
oscillatcr  cscillates  between  two  positions,  while  the 
whistler  must  oscillate  between  many  circumferential 
positions.  That  such  an  oscillation  takes  place  seems 
indicated  by  the  nozzle  configuration  and  the  separation 
reattach-Tient  cycle  noted  by  Hill  and  Greene. 

Seference  10  provides  a good  summary  of  the  mechanism  of 
the  Coanda  effect,  as  well  as  its  characteristics  regarding 
reattachment  and  sensitivity  to  sound.  Though  the  following 
comments  relate  specifically  to  two-  dimensional  flows,  it 
is  thought  that  the  general  behavior  and  mechanism  should 
carry  over  to  a three-  dimensional  situation. 


Figure  3 illustrates  the  Coanda  mechanism  in 
two-dimensional  flow.  As  shown,  a jet  is  issuing  from  an 
orifice  and  flowing  between  two  adjacent  walls.  Because  the 
jet  will  at  some  time  be  nearer  to  one  wall  than  the  ether, 
and  because  the  jet  is  entraining  air,  a'  pressure  difference 
across  the  jet  will  develop  that  pulls  the  jet  to  the  closer 
wall.  This . attachment  forms  a region  of  circulating  fld.w. 
The  point  of  jet  attachment  is  determined  by  this  region,  in 
which  a balance  is  struck  between  flow  entrained  by  the  jet 
and  flow  reinjected  by  the  jet. 
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Anything  that  disrupts  this  balance  causes  the  attachment 
point  of  the  jet  to  shift. 

Observation  has  shown  that  the  attachment  point  in 
Coanda  effect  is  sensitive  to  jet  velocity  and  sound,  within 
a range  of  Reynolds  numbers.  In  the  range  of  Reynolds 
number  from  200  to  3000,  under  the  effect  of  increasing  jet 
velocity,  or  a sound  input  of  increasing  amplitude,  the 
attachment  point  has  been  observed  to  move  upstream.  This 
upstream  movement  appears  to  approach  a limiting  position, 
beyond  which  the  point  will  not  move. 

Reference  10  relates  the  effects  of  jet  velocity  and 
sound  on  the  attachment  point  to  the  point  of  transition 
from  laminar  to  turbulent  flow  of  the  jet.  Both  increasing 
jet  velocity  and  a sound  input  will  move  the  transition 
point  closer  to  the  jet  exit.  This  movement  will  increase 
the  length  of  turbulent  flow  across  the  recirculation 
region.  A turbulent  jet  entrains  more  air  than  a laminar 
jet,  and  thus  the  balance  of  flow  in  the  recirculation 
region  is  disrupted,  with  a resulting  upstream  move  of  the 
attachment  point. 
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SCOPS  OF  GESEAKCH 


The  intent  of  this  research  'rfas  to  provide  a clearer 
insight  into  the  workings  of  the  whistler  nozzle,  and  to 
develop  an  improved  understanding  of  the  ways  in  which  one 
may  effect  control  over  the  characteristics  of  a free  jet. 


Specific  areas  of  research  were  the  nozzle  thrust 
efficiency,  the  nozzle  entrainment  rate,  and  visualization 
of  the  lip  interaction  process  described  by  Hill  and  Greene. 
Accordingly,  experiments  were  conducted  to  obtain  relevant 
data.  Thrusts  were  measured  on  a thrust  bed,  mass 
entrainment  measured  in  a device  patterned  after  that 
developed  by  Ricou  and  Spalding  (reference  7)  , and  high 
speed  motion  pictures  taken  of  tufts  about  the  nozzle  lip. 
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IV.  NOZZLE  DESCRIPTION 


The  whistler  nozzle  used  in  all  tests  is  illustrated  in 
figure  4.  It  consisted  of  a small  plenum/  convergent 
section,  short  constant  area  section  and  a moveable  collar 
over  the  constant  area  section.  The  inner  whistler  diameter 
was  1 inch,  while  the  collar  diameter  was  1.5  inches.  The 
nozzle  configuration  was  essentially  that  used  by  Hill  and 
Greene  in  their  investigations.  The  characteristic 
dimension  used  for  Reynolds  number  and  nondimensionalization 
of  axial  distance  from  the  nozzle  was  the  inner  nozzle 
diameter,  1 inch. 
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DESCRIPTION  OF  EXPERIMENTS 


A.  THRUST  EFFICIENCIES 


Thrust  efficiencies  were  obtained  by  measuring  nozzle 
thrust.  A comparison  was  made  between  measured  and 
calculated  thrust. 


Test  equipment  consisted  of  a thrust  bed  and  the  nczzle. 
The  thrust  bed  was  configured  to  measure  horizontal  thrust, 
and  slid  on  two  tracks,  mounted  to  a solid  foundation.  The 
nozzle  plenum  was  clamped  to  the  thrust  bed  which  was  in 
turn  secured  to  the  foundation  through  a lead  ring,  mounted 
with  a load  cell,  the  device  that  actually  measured  thrust. 
Air  was  supplied  to  the  nozzle  from  a large  tank,  connected 
to  the  plenum  through  a regulator  and  a flexible  tube  which 
led  vertically  away  from  the  plenum,  so  as  not  to  bias  the 
thrust  measurements  in  any  way.  Calibration  was  performed 
with  the  use  of  a tray,  connected  to  the  bed  through  a cable 
and  pulley  arrangement,  upon  which  one  placed  weights  of 
known  magnitude.  Peripheral  gear  required  included  a 
voltage  source,  digital  voltage  meter,  and  Wheatstone 
bridge.  Figure  5 illustrates  the  set’  up.  Plenum  pressures, 
were  measured  using  a mercury  manometer. 


Calibration  of  the  thrust  bed  indicated  a certain  amount 
of  hysteresis  in  the  system,  as  may  be  seen  in  figure  6. 
Not  serious  at  thrusts  of  about  6 pounds,  it  did  introduce 
problems  at  low  thrust  levels,  and  precluded  reliable  thrust 
measurement  at  low  Reynolds  numbers. 

Thrust  measurements  were  made,  for  various  collar 
positions,  over  a range  of  nozzle  exit  velocities.  Recorded 
data  included  gage  pressure  in  the  nozzle  plenum,  thrust 
readings  from  the  voltmeter,  collar  extension,  and  ambient 
pressure  and  temperature. 

Measured  thrusts  were  plotted  against  pressure  ratio 
across  the  nozzle.  Figure  7 presents  an  example  curves 
Ideal  thrust  was  calculated  according  to  the  procedures  in 
appendix  h.  Thrust  efficiencies  were  then  obtained  and 
plotted  against  nozzle  pressure  ratio,  exit  velocity,  and 
Reynolds  number  (based  on  the  nozzle  diameter) . 


B.  MASS  ENTEAIKflSNT 


Mass  entrainment  was  measured  directly,  using  a device 
patterned  after  one  developed  by  Ricou  and  Spalding 
(reference  7).  It  was  believed  that  this  device  should 
yield  more  accurate  results  than  an  approach  using 
integrated  velocity  profiles,  especially  in  the  case  of  an 
unsteady  nozzle.  Measurements  were  made  at  various  axial 
positions,  for  several  collar  positions. 


The  device  used  will  hereafter  be  referred  to  as  the 
entrainment  chamber.  It  was  designed  with  the  aid  of  data 
from  a research  report  by  Peschke  (reference  8)  . Appendix  B; 
provides  a detailed  account  of  the  design. 
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Plenum  Pressure/Atm.  Pressure 


The  principle  (f  operation  is  relatively  simple.  The 
entrained  air  of  a free  jet  flows  in  at  a right  angle  to  the 
jet  axis.  The  entrainment  chamber  simply  blocks  this  inflow 
with  a circumferential  membrane/  and  seals  the  ends  with  a 
plate  at  the  nozzle  and  an  orifice  at  the  axial  station 
where  one  wishes  to  measure  the  total  mass  flow.  Ideally/ 
all  of  the  entrained  air  is  supplied  to  the  exterior  of  the 
membrane  from  a regulated  source.  As  a free  jet  does  not 
support  any  axial  pressure  gradient,  the  proper  amount  of 
air  to  supply  through  the  membrane  is  indicated  by  the 
disappearance  of  any  pressure  gradient  across  the  orifice. 
One  may  regulate  the  membrane  flow  to  obtain  this  point. 
These  primary  aspects  of  the  device  are  illustrated  in 
figure  8. 

Critical  design  aspects  of  the  entrainment  chamber, 
especially  the  orifice  sizes,  are  treated  in  appendix  B.  It 
is  sufficient  here  to  simply  state  that  one  entrainment 
chamber  was  designed  to  perform  measurement  of  the  total  jet 
mass  flow  at  various  axial  locations. 

One  primary  difference  between  this  chamber’s  gecmetr-y^ 
and  that  of  ethers  (references  7 and  8)  should  be  noted.  As 
shown  in  figure  8,  exit  orifices  were  positioned  in  the 
porous  cylinder  at  various  distances  from  the  nozzle.  That 
portion  cf  the  porous  cylinder  extending  above  the  orifice 
plate  was  sealed  by  a fiberglass  cylinder  attached  to  the 
orifice.  Thus, in  this  situation  that  portion  of  the  jet 
beyond  the  orifice  and  inside  the  porous  cylinder  exit  was 
not  in  a truly  free-jer  flow  situation.  This  was  as  opposed 
to  the  reference  7 and  8 configurations  where  the  exit 
orifice  corresponded  to  the  exit  of  the  porous  cylinder. 
Thus,  in  those  experiments  the  jet  past  the  orifice  was  in 
truly  free-jet  flow  situation.  The  conf iguratiof  used  in 

this  report  was  adopted  in  order  to  facilitate,  the 
measurement  of  mass  flows  at  many  axial  locations.  It  was 


hoped  that  this  chamber's  geometry  would  not  affect  the 
chamber  readings  greatly,  although  some  effect  was  to  be 
expected . 

Peripheral  equipment  to  the  basic  entrainment  chamber 
included  flow  regulating,  metering,  and  measuring  devices. 
Figure  9 presents  a photograph  of  the  chamber  and  its 
peripheral  gear,  separate  air  sources,  both  fed  by  a large 
tank,  were  supplied  for  the  nozzle  and  membrane.  Nozzle 
flow  was  metered  through  an  orifice  and  controlled  with  a 
regulator.  Eecause  of  the  much  greater  mass  flow  reguired, 
membrane  air  was. supplied  through  a venturi,  followed  again 
by  a flow  regulator.  Pressures  across  the  orifice  were 
measured  with  an  inclined  water  manometer.  Pressures  across 
the  venturi  were  measured  with  a mercury  manometer. 
Pressure  across  the  entrainment  chamber  orifices  was 
measured  with  an  alchohol  raicromanometer . Nozzle  plenum 
pressure  was  measured  with  a water  manometer. 
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Membrane  air  temperature  was  measured  just  after 
venturi,  usin.g  a mercury  thermometer. 
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Entrainment  measurements  were  made  at  various  axial 
locations  for  various  collar  positions.  Eecorded  data 
included  pressure  difference  across  the  venturi  and  orifice, 
pressure  difference  across  the  entrainment  chamber  orifice, 
nozzle  plenum  pressure,  ambient  pressure,  air  temperature  in 
the  membrane  supply,  and  collar  position.  Standard 
procedure  for  each  run  was  to  vary  the  membrane  airflow  so 
that  one  obtained  both  positive  and  negative  pressure 
gradients  at  the  chamber  orifice.  This  provided  a more 
precise  location  of  the  proper  entrained  flow  rate.  It  is 
of  interest  to  note  that,  due  to  hardware  limitations  and 
the  large  aacunt  of  air  entrained  by  the  jet  at  the  furthest 
axial  locations  measured,  nozzle  exit  velocities  were 
limited  to  around  100  feet  per  second  at  the  further  axial 
locations.  This  also  affected  the  number  of  useful 
locations  for  whistler  measurements,  as  a strong  whistler 
oscillation  required  a relatively  large  exit  velocity,  say 
300  feet  per  second. 

Data  reduction  followea  standard  procedures  for  the 
venturi  and  orifice.  These  are  presented  in  appendix  C. 
Entrainment  was  plotted  as  M/MO {total  mass  flow  divided  by 
the  primary  mass  flow) , against  X/D,  axial  location  divided 
by  nozzle  diameter.  Working  curves  of  chamber  orifice 
pressure  gradient  are  presented  as  pressure  difference 
against  tt/liO,  in  the  example  plotted  in  figure  10. 


igure  10  - ORIFICE  PRESSDRE  DIFFERENCE  VS  EHTRAINMFNl? 
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C.  FLOW  VISUALIZATION 
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Flow  visualization  of  the  nozzle  jet  lip  interaction  was 
performed  in  two  ways:  high-speed  motion  photography  of  a 
ring  of  tufts  at  the  nozzle  lip/  for  the  purpose  of 
discerning  a pattern  to  the  interactions;  and  a mixture  of 
lamp  black  and  kerosene  inside  the  nozzle  collar  in  order  to 
observe  the  behavior  of  the  separation  reattachment 
location. 

Photcgraphic  equipment  consisted  of  a Honeywell  Eentax 
single  reflex  camera,  a Bolex  supreme  handheld  16rom  movie 
camera,  and  a Redlake  Hycam  16mm  high-.speed  movie  camera. 
The  tuft  ring  was  simply  constructed  of  cotton  threads  held 
by  a masking  tape  ring,  and  a photo  of  the  ring  is  presented 
in  figure  11.  Observation  in  the  lampblack  experiment  was 
performed  with  rhe  naked  eye  and  no  peripheral  equipment  was 
necessary. 

Photographic  tests  were  to  be  done  in  three  phases. 
High-speed  still  shots  were  first  taken  to  determine  the 
usefulness  of  the  tufts  for  visualization,  and  to  obtain  the 
approximate  required  film  speed  for  motion  photography. 
Next,  medium  speed  motion  pictures,  64  frames  per  second, 
were  taken  to  confirm  the  necessity  of  the  high-speed 
camera.  Finally,  high-speed  motion  pictures  were  taken  to 
visually  control  • the  tuft  motion  such  that  useful 
observations  could  be  made.  No  particular  collar  settings 
or  exit  velocities  were  used  in  photographic  work.  The  only 
requirement  was  to  produce  the  oscillation  phenomenon,  which 
was  never  difficult.  As  it  turned  out,  a lack  of  time 
prevented  the  completion  of  the  high-speed  photography 
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Lampblack  tests  were  performed  in  two  ways.  The  inside 
of  the  cellar  extension  was  liberally  coated  with  the 
lampblack  mixture,  and  either  exit  velocity  varied  at 
constant  collar  setting,  or  collar  setting  varied  at 
constant  exit  velocity.  Observations  were  all  visual,  none 
quantitative,  save  noting  the  range  of  exit  velocities 
covered,  which  was  from  about  100  to  300  feet  per  second.. 
The  lampblack  mixture  was  observed  to  pool  as  illustrated  inv 
figure  12,  and  the  downstream  — most  edge  of  the  pool  was 
interpreted  as  indicating  the  location  on  the  collar  of 
reattachraent  and  separation. 
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VI.  RESULTS  AND  DISCUSSION 


Nozzle  thrust  efficiency  results  are  presented  in  figure 
13,  plotted  for  various  collar  positions.  The  whistler 
oscillation  was  produced  at  each  collar  setting  for  the 
whole  range  of  Reynolds  numbers  shown,  although  sudden 
frequency  jumps  did  occur  in  this  range.  These  frequency 
jumps  did  not  appear  to  change  the  nozzle  thrust  or  thrust 
efficiency  in  any  discontinuous  manner.  Instead  the  trend 
was  as  indicated,  progressively  .decreasing  nozzle  efficiency 
with  increasing  collar  extension.  The  zero  collar  extension 
efficiencies  were  found  to  be  about  6 percent  lower  than 
those  found  in  a previous  report  (reference  8) . This  is 

believed  due  to  the  extremely  simple,  unflared  configuration 
of  the  nozzle,  which  may  be  seen  in  figure  4.  Reference  8 
nozzles  were  divergent  nozzles  and  as  such  should  be  more 
efficient.  Important  to  note  is  that  the  whistler 

oscillations  may  be  obtained  without  much  reduction  in  the 
basic  nozzle  efficiency. 

Entrainment  results  are  presented  in  figure  14,  plotted 
together  with  the  results  of  Hill  and  Greene  (reference  1) 
for  comparison.  As  may  be  seen,  in  the  case  of  the  basic 
jet  agreement  is  good  qualitatively  but  poor 
quantitatively.  In  the  case  of  the  whistler  nozzle 
agreement  is  not  good  in  any  way.  It  is  believed  that  the 
poor  guantitative  comparison  was  due  to  interference  of  the 
deep  orifice  cylinders  with  the  jet.  These  create  a 
pressure  gradient  from  the  orifice  to  the  chamber  exit  which 
quite  possibly  affects  the  readings.  As  noted  in  the 
experiment  section,  and  in  Appendix  C,  these  orifice 
cylinders  constitute  a major  difference  in  configuration 
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from  the  chambers  of  reference  7 and  8,  Additionally,  in 
the  case  cf  the  whistler  nozzle,  it  is  believed  the  method 
of  enclosing  the  jet  up  to  the  axial  position  of  measurement 
may  interfere  with  the  sound  - jet  interaction  (especially 
the  cloth  porous  cylinder)  and  the  formation  of  a lafge- 
scale  turbulent  structure.  Figure  15  illustrates  that 
readings  taken  were  relatively  constant  with  Reynolds 
number.  As  noted  by  Ricou  and  Spalding,  the  entrainment 
ratio  can  be  a strong  function  of  Reynolds  number  in  certain 
ranges  of  Reynolds  number-.  It  is  interesting  to  note  that 
this  behavior  is  supported  by  the  present  notion  of  jet 
turbulence  structure  and  its  development  as  exit  velocity 
increases,  which  postulates  that  the  major  contributor  to 
entrainment,  the  large  scale  structure,  dissappears  as 
Reynolds  number  increases. 


Figure  14 


MTRAIMHERT  RATIO  RESULTS  AND  OOMPAFiSON  DAM 


Lack  of  time  precluded  testing  at  more  jet  axial 
locations,  as  each  location  re-;aired  a different  orifice 
cylinder  (see  appendix  A)  . Of  the  original  four  orifice 
cylinders  constructed  using  reference  8 data,  only  two  were 
useable  as  the  spreading  angles  encountered  averaged  about 
28  degrees,  or  10  degrees  less' than  those  used  in  reference 
8,  In  order  to  be  effective,  the  orifices  had  to  be 
positioned  at  axial  distances  further  than  they  were 
designed  fcr.  One  of  those  was  modified  to  obtain  readings 
at  its  design  position,  also.  The  larger  two  were  never 
effective.  Another  problem  encountered  was  the  mass  flow 
limitaticn  cf  the  porous  cylinder  air  supply.  This 
restricted  the  useable  nozzle  flow  at  axial  locations  beyond 
X/D  of  about  6 to  such  low  values  that  a whistler 
oscillation  strong  enough  to  affect  the  entrainment  could 
not  be  produced. 

It  was  hoped  that  high  speed  motion  pictures  would  be 
taken  of  the  tuft  experiment.;  however,  this  phase  of  the 
photographic  experiments  was  not  reached.  T.ie  effectiveness 
of  the  tufts  in  indicating  the  circumferential  collar 
positions  of  attached  and  detached  flow  was  proven  in  tne 
high-speed  still  picture  phase.  Motion  pictures  with  a 64 
frame  per  second  camera  showed  that  a much  higher  speed 
camera  was  necessary  to  slow  the  motion  such  that  motion 
pictures  would  be  useful  in  slowing  the  tuft  motion  to  a 
reasonable  speed.  The  necessary  film  speed  was  estimated  as 
about  800  frames  per  second.;  however,  the  pictures  were  never 
taken.  Seme  representative  still  pictures  are  presented  in 
figure  16. 


Lampblack  tests  appeared  to  be  productive  as  indicators 
of  the  general  behavior  of  the  coll'ar  attachment  point. 
Indications  are  that  a sort  of  Coanda  effect  is  involved, 
modified  by  the  fact  that  while  the  jet  is  attached  at  any 
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point  on  the  collar  there  are  other  points  where  detached 
flow  allows  a backflow  into  the  collar  jet  gap,  which  tends 
to  relieve  the  Coanda  effect  causing  a circumferential 
shift  in  the  point  (or-  points)  of  attachment.  At  a fixed 
collar  extension,  with  increasing  jet  velocity,  the  attachment 
point  was  otserved  to  move  upstream  to  a limiting  position, 
just  as  the  Coanda  jet  described  in  reference  10.  However, 
at  fixed  exit  velocity  and  increasing  collar  extension  a 
different  behavior  was  observed.  The  attachment  point 
seemed  to  move  wrth  the  collar  until  a point  was  reached  at 
which  the  frequency  of  oscillation  jumped.  At  that  instant 
the  point  of  attachment  also  appeared  to  jump  suddenly  back 
to  its  original  position  relative  to  the  jet  exit.  This 
behavior  indicared  some  interdependence  of  the  whistler 
attachment  phenomenon  and  the  oscillation  phenomenon. 


In  addition  to  the  above  experiments,  two  other 
observations  were  made  which  verified  those  of  Hill  and 
Greene. 


First,  it  was  noted  that  extension  of  the  whistler 
collar  seemed  to  increase  mass  flov/  through  the  nozzle. 
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TUFT  PICTURES,  RESULTS 


This  was  confirmed  with  static  pressure  measurements  in  the 
nozzle  ’throat’.  These  indicated  that  the  collar  was  acting 
as  an  overexpanded  diffuser,  resulting  in  an  effective  exit 
area  slightly  greater  than  the  ’throat’  cross  section  but 
less  than  the  collar  cross  section,  and  further  resulting  in 
a greater  mass  flow  through  the  nozzle. 

Second,  it  was  noted  that  a correspondence  between 
freguency  of  oscillation,  jet  exit  velocity,  and  collar 
extension  seemed  to  exist.  It  was  found  that  by  gradually 
extending  the  collar  as  plenum  pressure  was  increased,  one 
could  ’follow’  a frequency  of  oscillation.  . Accordingly, 
collar  extensions  were  measured  and  jet  velocity  calculated 
at  fixed  frequencies  of  oscillation,  and  it  was  found  that 
these  fixed  frequencies  corresponded  to  constant  ratios  of 
collar  extension  to  jet  velocity,  thus  confirming  the  organ 
pipe  character  of  the  whistler  oscillation  noted  by  Kill  and 
Greene. 


VII.  CONCLUSIONS  AND  RSCOMHENDATIONS 


The  thru53t  efficiency  and  mass  entrainment 
characteristics  of  the  whistler  nozzle  have  been  measured, 
in  the  case  of  the  entrainment,  using  a technique  heretofore 
not  used  on  an  unsteady  jet.  Further,  observations  have 
been  made  which  should  contribute  to  a further  understanding 
of  the  whistler  flow  mechanism. 

Several  directions  of  further  work  are  indicated. 
First,  it  is  recommended  that  a thrust  measurement  device 
using  a pendulum  arrangement,  similar  to  that  used  in 
reference  15,  be  constructed  for  a more  accurate 
determination  of  the  whistler  efficiency.  A device  of  this 
type  might  avoid  the  hysteresis  problems  encountered  with 
low  thrusts  on  a conventional  thrust  bed.  It  is  also 
recommended  that  a whistler  nozzle  of  smaller  throat 
diameter  be  constructed  to  permit  mass  entrainment 
measurements  , at  further  jet  axial  locations,  in  the 
entrainment  chamber,  at  the  high  exit  velocities  required  to 
produce  a strong  oscillation.  This  will  also  permit  further 
study  of  the  possibility  of  whistler  measurements  in  the 
chamber  ever  a greater  range  of  Reynolds  number.  Next  it  is 
recommended  that  further  flow  visualization  experiments  be 
conducted.  Specifically,  high  speed  motion  pictures  of  the 
tufts  used  here,  and  high  speed  schlieren  photography  of  the 
jet,  should  prove  useful.  Finally,  it  is  recommended  that 
other  unsteady  nozzles  be  tested  in  the  chamber,  as  the  poor 
results  obtained  here  with  the  whistler  are  only  preliminary 
and  may  not  be  indicative  of  rhe  usefulness  of  the  chamber 
with  non-acoustic  effects  such  as  swirl. 


' VIII.-  INIROPUCTION-  THE  RING  WING 
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I The  present  work  concerns  a hollovr  cylinder,  axxally 

I aligned,  in  supersonic  flow  and  undergoing  small  amplitude 

■j  oscillations  in  angle  of  attack.  The  intent  is  to  calculate 

■j  the  resulting  pressure  distributions  and  aerodynamic  mcnieuts 

on  what  might  be  more  properly  termed  a ring  wing. 


Previous  investigations,  i.e.  references  12  and  13,  have 
applied  the  method  of  characteristics  towards  the  separate 
calculation  of  the  inner  and  outer  flowfields,  as  well  as 
pressure  distributions  and  generalized  aerodynamic  forces, 
for  such  a cylinder  in  supersonic  flow  whose  walls  are 
undergoing  small  amplitude  sinusoidal  oscillations,  or  panel 
flutter.  These  previous  investigations  resulted  in  computer 
programs  to  perform  the  necessary  calculations. 


The  primary  thrust  of  the  present  work  was  the 
modification  of  the  previously  created  programs  and  then  the 
ccmbinaticn  of  their  results  to  calculate  the  unsteady 
pressure  distribution,  lift  and  moment  on  the  cylinder 
undergoing  angle  cf  attack  oscillations.  Finally,  using  the 
modified  programs,  a parametric  study  of  these  items  was 
conducted. 
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IX.  PROBLEM  FORMULATION 


This  development  is  closely  based  on  several  previous 
papers  , i.e.  references  12  and  13. 

Consider  a circular  cylindrical  shell,  in  supersonic, 
inviscid,  adiabatic  £1oh,  whose  axis  is  aligned  with  the 
freestreara,  and  which  is  undergoing  small  amplitude 
oscillations  in  angle  of  attack,  about  the  Z axis,  (figure 
17) 

The  governing  equation  for  this  flow  situation  is  the 
linearized  unsteady  potential  equation,  in  cylindrical 
coordinates. 


(1-mM  «ee  - 2M  «xt  ■ h *'tt  = ° 

3T  3T  O C 


The  relevant  boundary  conditions  are  the  flow  tangency 
condition  and  Sommerfeld* s radiation  condition,  i.e., that 
disturbances  will  propagate  away  from  their  source. 


Jvlpha,  Oscilation  Amplitude 


ss 

fete 


tte:'^ 

If^'- 

ll  :t  ,.i 

i •: 

|i : f:  v" 

P •!%  '; 

!:#■ 
V -'-Si..-, 


One  further  condition  is  necessary^  along  the  X-^axis,  which 
■aiaounts  to  setting  the  pressure  disturbance  at  the  axis 
equal  to  zero. 

The  linearized  form  of  the  flow  tangency  condition  is: 

^rl  = Jji  + U 9h  ,0  ■$:  X ;$  L 
r=R  3t  <»  9x 


which  must  be  satisfied  on  the  cylinder  inner  and  outer 
surface. 

As  stated,  the  boundary  condition  at  the  axis  is  a zero 
pressure  disturbance: 


CpI  = 0 


The  linearized  form  of  the  pressure  coefficient  is  as 
fellows: 


Cp  = -2  - A.  = P - Pc 

^ ^ ^ Uco  q„ 


After  nondinensionalization,  using  cylinder  length  and 
length  divided  by  freestream  velocity,  the  above  equations 
become: 


' -i' 
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2M^  ♦tt  = ° 


Cp  = -2(  *j.  + 4.^  ) 


subject  to: 


*J’rl  “ — /O  ^ 1, 


r=R  3t  3x 


One  now  assumes  a cylinaer  oscillation  of  the  form: 


h(x,e,t)  - Z(x)  cos  (8)  e 


ikt 


where  the  nondiraensional  frequency,  k,  is; 


k = 0)  1 

U„, 


The  deflection  amplitude,  Z (x) , for  the  case  under 
consideration,  is  of  the  form: 


Z (x)  = a(  X - X ) 
0 
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is 


Following  the  assumed  form  of  the  deflection,  the 
perturbation  potential  must  be  of  the  form: 


(Ji(x,r,e,t)  = <|)(x,r)  cos(6)  e 


ikt 


and,  after  substitution,  the  basic  equations  become: 


Q-mM  d>....  + (b  + 1 (!)  + (k^M^-l  ) (j)  - i2kM  <j).,  = 0 

-*'xx  ■»'j-r  — r zrr-  ^ 


Cp  = -2(  (f)j^  + ik<i>  ) 


subject  to: 


(J)  1 = Z + ikZ  ,0  $ X 1. 

^ «PS  ^ 

r=R 


and; 


C = 0 

P r=0 
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X.  «3TK0D  OF  CHARACTERISTICS 


For  supersonic  flow,  the  linearized  perturtation 
potential  equation  is  hyperbolic,  and  has  characteristics 
which  satisfy  the  following  differential  equation: 


(M^-l)  d^r  - d^x  = 0 


Since  the  characteristics  in  supersonic  flow  correspond 
to  the  Mach  lines,  then  for  ds,  the  differential  arc  length 
along  a characteristic; 


^ = + 1_ 
ds  M 


dx  = /M^-1 
ds  M 


An  arbitrary  function,  F(x,r),  has  the  following 
derivatives  along  the  characteristics: 


dF  = F„  dx  + F^  dr  = /m^^I  F^  + 1 F 

as  ^ as  ^ di  n~  M ^ 


Defining  dsl  and  ds2  to  be  the  differential  arc  lengths 
along  the  left, and  right-running  Mach  lines  respectively, 
then; 


Fi  = 


, F,  = dF 


1.45^  ' 'J 
f.  '%■.  ' 


4. ’a.  v,''2 


, t_  S'.  ~Z.  *i.r  '' 

fVt 


t >,  '-I  .;<*5 


fjjf 


. dSi 


The  cross  derivatives  become: 


-F^-  = = 1 (d-M^-)  P + P ) 

12  21  ' XX  rr 


Solving  now  for  Fx  and  Fr; 


F^  = (F.  + F^) 

^ 2/M-l  ^ ^ 


Pr  = (P^  - Pj) 


Now  the  basic  equations  may  be  written  in  terms  of  <j) 
and  its  derivatives  along  the  characteristics  as  follows: 


>12  ""  1 ) <l>  " AISM-  (<1)^  + 

21  2rM  ^ ^ M^-1  1 2 


<j)j.  = ii  ((f)^  - <l>^)  = + ikZ 

2 


C = -2  (4)  + ik<i>) 

P X 


These  equations  have  been  put  into  finite  difference 
form  and  programmed  in  Fortran  IV  for  the  separate  cases  of 
exterior  flow  and  interior  flow  of  a cylinder  whose  walls 
were  undergoing  small  amplitude  sinusoidal  oscillations.  In 
the  interior  flow  program  the  X axis  boundary  condition  is 
satisfied  on  a small  inner  cylinder  in  place  of  the  axis^  to 
circumvent  problems  encountered  should  r go  to  zero.  For 
detailed  acccunts  of  the  equations  and  programs  refer  to 
references  12  and  13. 


The  primary  work  here  was  the  transformation  of  zhe 
aforementioned  programs  to  the  case  under  consideration  , 
i.e..  a cylinder  undergoing  small  amplitude  oscillations  in 
angle  of  attack,  and  then  a parametric  variation  study  of 
the  pressure  distributions  and  generalized  aerodynamic 
force.  The  programs  were  modified  such  that  all  cases  run 
were  with  an  angle  of  attack  amplitude  of  0.1  radians. 
Since  the  basic  equation  is  in  a linearized  form,  the 
conversion  of  any  results  to  other  angle  of  attack 
amplitudes  is  a simple  matter  of  multiplication.  For  those 
interested,  the  programs  resulting  from  changes  required  to 
the  reference  12  and  13  programs  are  listed  in  appendix  D. 


XI.  grnemlized  ME22IM11IS  force 


As  stated  above,  the  results  studied  were  the  pressure 
coefficient  amplitude  distribution,  and  the  generalized 
aerodynamic  force,  Qmr.  In  its  most  general  form,  Qmr  is  as 
follows: 


Qmr  = 1 /C  (x)  dx 

y Pm  ^ 

Z(x)  = Z A.  'i'.  (x) 

j ^ ^ 


where  Cpm  is  the  pressure  coefficient  amplitude  resulting 
from  the  ra’th  axial  deflection  mode,  PSIr  is  the  r'th  axial 
deflection  mode,  and  A is  the  amplitude  of  the  j’th  axial 
deflection  mode.  For  the  special  case  Under  consideration, 
Q takes  the  following  form: 

Q = 1 ;c  'l'{x)  dx 

T P 


Z (x)  = (x)  = a(x  - x^) 


Thus  , in  the  case  under  study  q has  the  special 
rpretaticn  of  being  the  ■ moment  coefficient  amplitude 
t the  point  Xo  multiplied  by  the  angle  of  attack 

itude,  alpha,  which  is  egual  to  0.1  radians  in  all  cases 
ented  . 


XII.  RESULTS  AND  DISCUSSION 


Several  cases  cf  comparison  were  considered  towards  the 
evaluation  of  program  results.  Figure  18  presents  two  cases 
of  comparison  for  pressure  distributions.  The  outer  radius'- 
to-  length  ratio  was  taken  to  such  a value  (10)  as  to 
approximate  a flat  plate  locally.  The  inner 
radius-to-length  ratio  was  taken  to  such  a value  (9.75)  as 
to  cause  a reflection  on  the  inner  surface  of  the  outer 
cylinder. 

For  sufficiently  large  radius  - to  - length  ratios  the 
exterior  pressure  distribution  should  approximate  the 
theoretical  value  for  a flat  plate  in  supersonic  flow,  i.e. 
2a/ ^ Figure  18-A  presents  this  comparison,  and  as  can  be 
seen  agreement  is  good. 


For  sufficiently  large  radius  - to  — length  ratios  the 
interior  pressure  distribution  should  approximate  the 
theoretical  value  for  a flat  plate.  Additionally,  if  the 
inner  cylinder  is  of  such  a radius-to-length  ratio  as  to 
cause  a Nach  wave  to  reflect  upon  the  outer  cylinder  inner 
surface,  a flow  situation  exists  which  should  approximate  a 
flat  plate  in  a free  jet,  where  the  mach  waves  are  reflected 
from  the  jet  boundary.  This  is  due  to  the  boundary 
condition  imposed  on  the  inner  cylinder,  Cp  = 0.  Reference 
16  has  treated  this  case  in  analytic  . form  for  small 
freguencies.  Figure  18,  A and  B,  presents  a comparison  of 
program  results  and  the  results  of  reference  16.  As  can  be 
seen  agreement  is  good. 


Reference  14  has  treated  the  case  of  a nonoscillating 
ring  wing  in  supersonic  flow.  Figure  19  presents  a 
comparison  of  program  results  and  those  of  reference  14,  for 
the  case  shown.  As  can  be  seen  agreement  is  fair. 


Tn  reference  12,  which  developed  the  basic  inner  flow 
field  program  used  in  this  report,  the  occurrence  of  a 
numerical  instability  in  the  pressure  distributions,  under 
certain  conditions,  was  noted.  It  was  found  that  this 
disturbance,  which  occurred  only  after  a reflection,  was  due 
to  an  undefined  constant  in  the  subroutine  which  computed 
the  boundary  condition  at  the  axis.  After  defining  the 
constant  correctlj'^  the  pressure  distributions  were  seen  to 
settle  down  quite  nicely.  Figure  20  presents  a before  and 
after  comparison  of  pressure  distributions  for  a 
representative  case. 


Reflection 


R/L  =.4 

rf./L=. 00001 

M =1.5 
k =0.0 


Figure  20 


N.UMEEICAL  INSTABILITY  ILLUSTRATION 


Figures  21  through  23  present  various  parametric  studies 
of  the  moment  coefficient,  Q,  for  cases  both  with  and 
without  a reflection  on  the  inner  cylinder  surface.  In 
these  figures  the  coefficient  plotted  is  actually  the 
difference  between  that  calculated  for  the  inner  cylinder 
surface  and  that  calculated  for  the  outer  cylinder  surface, 
divided  by  the  angle  of  attack  amplitude.  Thus  the 
coefficients  plotted  are  total  moment  and  pressure 
coefficients,  as  opposed  to  those  plotted  earlier,  which  are 
the  coefficients  for  the  inner  and  outer  cases  separately. 
As  can  be  seen  from  figure  22, for  a center  of  rotation  at 
the  leading  edge  the  cylinder  has  definite  geometric 
regions  of  stability  with  respect  to  angle  of  attack. 
Figure  21  illustrates  the  effect  of  increased  frequency  on 
the  moment  coefficient.  Figure  23  illustrates  a typical 
pressure  distribution  including  one  reflection,  front  which 
one  can  infer  the  reason  for  the  behavior  of  Q in  figure  22, 
and  also  the  behavior  of  Q with  a varying  center  of 
rotation. 
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Figure  21 


PITCHING  H0N2NT  VS  RSDUCED  FREQUENCY 


XIII 


CONCIOSIOWS  AND  RSCOHttENDATIO’.IS 


The  computer  programs  of  references  12  and  13  have  been 
successfully  adapted  to  the  case  of  a ring  wing  in 
supersonic  flow  undergoing  oscillations  in  angle  of  attack. 
Available  cases  for  comparison  have  shown  good  agreement, 
indicating  the  validity  of  the  method.  The  programs  make 
possible  the  relatively  easy  calculation  of  the  unsteady 
lifts  and  moments  of  a ring  wing  in  supersonic  flovr,  for 
such  applications  as  tubular  artillery  shells  {STOP, 
spinning  tubular  projectile)^  although  spin  may  induce 
aerodynamic  effects  not  accounted  for  in  the  method. 
Further  evaluation  and  parametric  studies  ace  recommended 
for  a -more  complete  verification  of  the  computer  programs. 


&PPEHDIX  A 


TSENTROPTC  THRUST  CALCULATIONS 


Xsentropic  thrust  calculations  made  use  of  standard 
isentropic  relations,  calculations  were  based  on  plenum 
pressure  and  nozzle  exit  area.  The  formula  for  isentroplc 
thrust  is  arrived  at  as  follows.  Nozzle  thrust  may  be 
written  in  the  following  form: 


= W V =(pAM^jrRT) 
ex 


ex 


where  M is  nozzle  mass  flow,  A is  exit 
applied  the  speed  of  sound  relation. 


area,  and  one  has 
Exit  HacH  number  may 


be  written  in  terms  of  pressure  ratio  as  follows: 


M 


ex 


= _2_  ((Pq/P  )^'Vy  -1) 


y-1 


ex 


Applying  this  and  the  state  equation  to  the  formula  fo: 
Ts  one  obtains: 


-1 


T = p 


ex  ex 


2 ((Pq/P  ) 


T=1 


ex 


'D 


Now,  once  one  has  measured  plenum  and  ambient  pressure, 
and  thrust,  Texp,  one  can  calculate  the  isentropic  thrust 
efficiency: 


.APPENDIX  B 

ENTRAINMENT  RATIO  CALCHLATIONS 


Mass  flows  and  entrainment  ratios  were  calculated 
according  to  the  equations  and  procedures  sat  out  in 
reference  11,  the  A5ME  Power  Test  Codes,  The  equation  for 
mass  flow,  from  reference  11  is; 


w = 359  c f Fa  Y /hp  ,lb/hr 


where; 


c Discharge  Coefficient 

^ =1//1-B‘*  e=d/D 

d and  D d.'^  | 

^ ) i 

Ea  Thermal  Expansion  Factor 

y (Ya  for  .1  venturil  Adiabatic  Expansion  Factor 
h The  pressure  drop  across  the  device,  in  inches  of 


water 

p Density  of  the  fluid  being  metered,  pounds  per 
cubic  foot 


If  w1  is  the  entrained  flow  and  w2  the  primary,  then  (w1 
uses  a venturi,  w2  an  orifice) : 

W-]  = c,  f,  (d, ) ^ Fa-,  Ya^  /h-, 

«2  T-2^  Ti^ 
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Non-geometric  factors  for  the  primary  flow,  such  as 
discharge  coefficient,  could  be  reasonably  approximated  for 
all  tests  as: 


C2=.65 


Fa2=1,0  Y2=.99 


Tn  the  case  of  the  entrained  flow  these  factors  had  to 
be  changed  for  each  test  axial  location,  due  to  the  widely 
varying  flow  requirements.  The  geometric  factors  above  had 
the  followino  fixed  values: 


F1=1.0328  d1=.75 


F?.=  1.0792  d2  = 1.11 


How,  entrainment  ratio  is  commonly  written  as: 

_M  = w.  + w = w,  + 1 

Mq  -W2 

Using  the  previous  formula  for  u1/w2,  the  entrainment 

ratio  may  now  be  written  as: 

M = A +1 

Mn  Tiu 


where  A is  the  combination  of  the  ratios  of  factors 
previously  listed.  This  aquation  was  used  to  calculate  the 
nozzle  entrainment  ratios. 


APPENDIX  C 

ENTRATNK3KT  CHAHBEP  DESIGN 


The  entrainment  chamber  design  was  closely  based  on 
previous  work  by  Ricou  and  Spalding  {reference  7) , and 
Peschke  (reference  8).  Peschke’s  report  especially  was  a 
good  source  of  specific  design  data,  while  Ricou  and 
Spalding  was  primarily  a source  of  the  design  principles. 

The  entrainment  chamber  consisted  essentially  of  a 
central  porous  cylinder,  closed  at  one  end  with  a solid  wall 
through  which  the  nozzle  protruded,  the  other  end  being  open 
to  the  outside.  Orifice  plates  of  specific  diameters  were 
constructed  and  mated  to  cylinders  of  specific  length, 
length  and  diameter  specified  according  to  the  data 
available  in  Peschke.  These  orifice  cylinders  were  of  such 
outside  diameter  that  they  fit  closely  inside  the  porous 
cylinder,  into  which  they  were  inserted  for  the  measurement 
of  mass  flow  at  the  various  jet  axial  positions 
corresponding  to  the  orifice  diameters.  The  purpose  of  the 
orifice  was  to  seal  the  open  end  such  that  entrained  air 
could  only  enter  through  the  porous  cylinder.  The  purpose 
of  the  orifice  cylinder  was  to  seal  off  that  portion  of  the 
porous  cylinder  beyond  the  jet  axial  position  where  flow  was 
being  measured.  The  porous  cylinder  was  surrounded 
everywhere,  except  at  the  open  end,  by  a box,  also  encasing 
the  nozzle  plenum,  which  served  as  the  membrane  air  supply. 

The  entrainment  chamber  was  about  4 by  4 feet  square  and 
4 feet  high,  supported  on  3 foot  legs.  The  porous  cylinder 
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was  30  inches  in  diameter  and  25  inches  long.  Orifice 
diameters  used,  corresponding  to  the  listed  jet  axial 
positions,  were  as  follows: 

Diameter  X/D 


1.75 


® .-M-- 

1-4=^ 


wsr  0>') 


y 


13.0 


The  entrainment  chamber  box,  legs,  orifice  plates  and 
porous  cylinder  frame  work  were  all  constructed  of  wood. 
The  orifice  cylinders  were  made  of  fiberglass  sheet  and  the 
porous  covering  to  the  porous  cylinder  frame  work  was  two 
layers  of  dacron  crepe  cloth,  as  suggested  in  Peschke’ s 
report. 

One  critical  chamber  design  parameter  was  orifice 
diameter.  It  is  clear  that  too  small  an  orifice  would 
interfere  with  the  jet  being  measured,  while  too  large  an 
orifice  would  not  properly  restrict  the  entrained  flow  to 
entering  only  through  the  membrane,  thus  preventing  the 
creation  of  a readable  pressure  gradient  across  the  orifice. 
The  approach  to  this  problem  in  both  references  7 and  8 was 
to  find  an  optimal  orifice  size  at  a specific  jet  axial 
location  by  testing  several  orifice  sizes.  To  minimize 
orifice  construction,  this  investigation  took  a different 
approach.  That  was  to  use  the  reference  8 orifice  data  as  a 
starting  point,  and  vary  the  orifice  axial  location,  such 
that  a readable  pressure  gradient  was  created  while  not 
interfering  with  the  jet.  in  both  approaches  the  attempi  as 
to  optimize  the  orifice  diameter,  jet  axial  location 
combination.  ' 


One  primary  difference  between  this  chamber's  geometry 
and  that  of  others  {references  7 and  8)  should  he  noted.  As 
shown  in  figure  8,  exit  orifices  were  positioned  in  the 
porous  cylinder  at  various  distances  from  the  nozzle.  That 
portion  of  the  porous  cylinder  extending  above  the  orifice 
plate  was  sealed  by  a fiberglass  cylinder  attached  to  the 
orifice.  Thus  in  this  situation,  that  portion  of  the  iet 
beyond  the  orifice,  and  inside  the  porous  cylinder  exit,  was 
not  in  a truly  free  jet  flow  situation.  This  was  as  opposed 
to  the  reference  7 and  8 configurations  where  the-  exit 
orifice  corresponded  to  the  exit  of  the  porous  cylinder. 
Thus  in  those  experiments  the  jet  past  the  orifice  was  in  a 
truly  free  jet  flow  situation.  The  configuration  usee)  in 
this  report  was  adopted  in  order  to  facilitate  the 
measurement  of  mass  flows  at  many  axial  locations.  Tt  was 
hoped  that  this  chamber's  geometry  would  not  affect  The 
chamber  readings  greatly,  although  some  effect  was  to  be 
expected. 

An  i.iaportant  problem,  again  with  regard  to  the  chamber 
orifices,  concerned  the  measurement  of  the  pressure 
differential  across  the  orifice.  Both  reference  7 and  8 
indicated  that  a measurement  accuracy  of  .001  inches  of 
water  was  required.  Reference  7 used  a micro  manometer 
while  reference  8 used  a pressure  transducer  for  this 
purpose.  After  the  entrainment  chamber  was  constructed 
experience  showed  that  such  accuracies  were  required.  This 
investigation  used  an  alchohol  micro  manometer  which 
provided  an  accuracy  of  .0001  inches  of  manometric  fluid. 
This  accuracy  was  obtained  through  the  combination  of  an 
inclined  tube,  hairline  indicator,  magnifying  lens  and  a 
micrometer  adjustment  to  lndicat.or  level. 

Flow  measurement  devices  for  the  primary  and  entrained 
flow  were  an  orifice  and  venturi  respectively.  These  were 
constructed  and  used  according  to  ASM3'  Power  Test  Codes 
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specifications  (reterence  ‘j'!). 
specifications,  an  inclined  water  manoiseter  was  found  to  be 
necessary  for  measuring  the  orifice  pressure  drop,  while  a 
vertical  mercury  manemetec  was  sufficient  for  the  venturi- 
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RPPBKDIX  D 


PROGRAM  LISriNGF,  AND  DOCUMENTATION 


This  appendix  contains  basic  program  and  input 
descriptions  for  the  two  programs  used,  from  references  12 
and  13.  The  programs  also  are  listed  in  their  modified  form 
for  the  computation  of  a ring  wing  at  angle  of  attack.  The 
program  statements  relevant  to  the  modification  are 
indicated  in  the  program  listings.  Also  indicated  in  the 
interior  flowfield  program  is  the  statement  whose  omission 
caused  the  numerical  instability  noted  in  reference  12. 

The  inner  flow  field  program  consists  of  8 subroutines, 
and  a main  program  which  manipulates  each  of  these.  The 
action  of  each  particular  subroutine  is  indicated  in  the 
program  listing.  Inputs  are  made  by  way  of  a Fortran 
Namelist.  The  input  parameters  are,  in  proper  order: 

FSTRIIN  Free  stream  mach  number. 

REDFRQ  The  reduced  frequency. 

RO  The  outer  cylinder  radius  to  length  ratio. 

El  The  inner  cylinder  radius  to  length  tio. 

NGRDFN  Grid  fineness,  or  the  number  of  grid  points 
taken  on  the  first  right  running  mach  line. 

K Circumferential  mode  number. 

H Axial  mode  number. 

Ml  The  m in  Qmn,  the  generalized  aerodynamic  force. 


7i| 


IPRTNT 


IPRTNT  An  output  parameter.  0 indicates  only  basic 
output.  1 indicates  output  including  such  things  as  Phi  and 
its  derivatives  at  each  flow  field  point. 


Of  the  above  inputs,  in  the  modified  program,  the  inputs 
N,  M,  and  HI  must  be  set  to  1,  as  the  case  under 
consideration  is  not  panel  flutter  but  a cylinder  at  angle 
of  attack.  The  other  inputs  are  arbitiary. 
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OOO  O OOO  O OOO  OOO  OOO  OOOO 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C C 

C CYLINDRICAL  SHELL/PANEL  FLUTTER  PROGRAM  C 

C C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C THIS  IS  THE  MAIN  PROGRAM 

C 

complex  PHn400J  ,DPHIS1(400)  ,0PHIS2(400) 

DIMENSION  X(400) » R (400 ) , XPT (400 ) ,CPR( 400 ), CPI (400) 

0 1 ME  NS  I GN  DUMMY ( 400 ) t OUNCE ( 400 ) 

CCMM0N/8LCK1/  FSTRMN , DELTAS t DXt  DH, DSTST R , RO  ,R I , REDFRO 
CCMMON/ BLCK2/  A 1 » A2 , A3 , A4 , RN, RH, RMl , FACT  1 , F ACT2 
CQMM0N/BLCK3/  NGROFN,NPTS, LCCUNT t IHAVEP  » N , M , Ml , I PR INT 
CCMM0N/BLCK4/  PH  I t DPHI SI , OPH IS2 
CGMM0N/BLCK5/  X, R,XPT  ,CPR,CPI 
DATA  PI/3.141593/ 

C 

ENDPTJY2,Y1,X0,X2,DELX)  = ( Y2  ~ Y1  )*(  XO  -X2  )/DELX 
1 + Y2 

C 

EPS  = l.E-04 
lER  = 0 

CALL  INPUT(IER) 

IF  ( lER  .EQ.  1 ) GO  TO  12 
1000  CONTINUE 
C 

I  HAVE P = 1 
LCCUNT  = 1 
C 

CALL  INITAL 
C 

1 CONTINUE 
C 

CALL  COMPXR 

IS  THE  GRID  POINT  ON  THE  INITIAL  RIGHT  RUNNING  MACH 
LINE? 

IF  ((  LCOUNT  .EQ.  1 ) .AND,  (IHAVEP  .LE.  NGRDFN  )) 

1 GO  TO  3 

IS  THE  GRID  POINT  ON  THE  OUTER  CYLINDER  SURFACE? 

IF  ( IHAVEP  .EQ.  0 ) GO  TO  5 

IS  THE  GRID  POINT  ON  THE  INNER  CYLINDER  SURFACE? 

IF  ( IHAVEP  .EQ.  NGRDFN  ) GO  TO  6 
THEN  THE  GRID  POINT  IS  IN  THE  GENERAL  FLOW  FIELD. 

2 CALL  GENFPT(IER) 

IF  ( lER  .EQ.  1 ) GO  TO  12 
IHAVEP  = IHAVEP  + 1 
GO  TO  1 

3 CALL  MACHLN 

IF  { IHAVEP  -.EQ.  NGRDFN  ) GO  TO  4 
IHAVEP  = IHAVEP  *■  1 
GO  TO  1 

THIS  IS  THE  LAST  POINT  ON  THt  INITIAL  MACH  LINE. 

4 LCCUNT  = LCOUNT  I 
IHAVEP  = 0 
6C  TO  1 

5 CALL  RAD2 

IS  THIS  GRID  POINT  PAST  THE  FLEXIBLE  CYLINDER  LENGTH? 
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ooou 


IF  ( {X(1)+EPS)  .GE. 
IHAVEP  = IHAVEP  *■  1 
GG  TO  1 


1.  ) GO  TO  8 


C 

C 


C 

C 

C- 

c 


c 

c 


CALL  RADI 
IHAVEP  = .0 
LCCUNT  = LCOUNT 
60  TO  1 


+ 1 


8 CONTINUE 

IF  ( ABSIXPTILCOUNT)  - 1.  ) ,LE.  EPS  ) GO  TO  10 
L = LCOUNT 
LMl  = LCOUNT  - 1 
LM2  = LCOUNT  - 2 
XFT<L)  = U 

CPR(Ll  = ENDPT(CPR(LMU  tCPR(LM2l  tl.OOtXPTILMU  tOSTSTR) 
CPKU  = ENDPT(CPI  (LMU  tCPI  (LM2)  ,1.00tXPT(LMl),0STSTR) 

WRITE(6,105» 

DO  9 I = ItL 
ARG  = RM1#PI*XPT( I ) 

THE  FOLLOWING  TWO  CARDS  WERE  MODIFIED  FROM  THE 
ORIGINAL  REFERENCE  12  PROGRAM. 

DUMMY!  I)=CPR{n*-l'-!=XPT(I) 

9 DUNCE! I )=CPIl I)*.1*XPT( I) 

SUMl  = 0.5*!CPR!L)  *■  CPP ! LHl ) )*!  1.00  - XPTILMD) 

SUM2  = 0.5T<(CPI!L}  i-  CPI  (LMl  ))#(  1.00  - XPHLMD) 

WRITE! 6t 1 14)  FSTRMNt  REDFRQ, RO, R I »N 

CALL  QSF ( OSTSTP , DUMMY , DUMMY , LMl ) 

CALL  QSF(DSTSTR, DUNCE»DUNCE,LM1) 

GREAL  = 0.5*(-DUMMY{LMU  + SUMl) 

QIMAG  = 0,5*(DUNCE(LM1)  + SUM2) 

WRITE!  6t  115)  MfMWGR6AL,M,Ml,>QIMAG 
GO  TO  12 

10  CONTINUE 
XP7f(L)  = X. 

WRXTE!.6t'10'5) 

DO  VX  I = WLCOUINT 
ARG  :=  RMl=''!^P-r-’=XPTin  )• 

DUMMY!  I,).=  CPR ! I ) ^^..J*XPT(  I ) 

1 1 DUNCEsC  I ) -CP  Id)*  .X*XPT ! I ) 

WRrTEr(v6.,  1X4)  FSTRMNt  PEDFRQ,  RO,  RI  tN 

CALL  QSjfaLOSTSTR,  DUMMY  tDUMMY  , L) 

C ALL  Q^Fl^DSrSTR  » OUNCE  t DUNCE  rL) 

QREAL  = 0i.'5*0UMMY(L) 

QXMfG=  0.-5*OUNCE(L) 

WRITE-.!6,tMt5.}  M,Ml,QReAL,M,Ml,QIMA6 

12  CONTINUE 
WRITEX6,99-g))i 

998  PORM'AT'!  • X'V^CP'  R E A L * ) 

CALL  PLOT PjCXtP-rrC.PR, LCOUNT , 0 ) 

.WP;XT&(=6,j99-'5f^ 

999  PORMA-X!  < X»V*''CP  I MAG  I NARY»  ) 

CALL  PLOTP  (:>CPr^ie P I , LC OUNT , 0 ) 

STOP 

150  FCPHA-T(T10.5.) 

105  FORMATdX*  rT19,»  I*  ,T27t'X*  ,T38,  »CPR»  ,T50  , • CPI » »// ) 

110  FGRMATdO*  tX5X,  I5,3!3X,F9.5  ) ) 

114  FORMWW/*  *,tT20,*MACH  NUMBER  = ‘tFlO.Tt//'  »»T20t 

1 'REDUCED  FREQUENCY  = ',F10.7,/‘  • ,T20 ,'! RADIUS /' , 

2 'L'ENCTH}  CUTER  = ',Fl0.7t/'  • »T20 1 '(  RADI  US/LENGTH ) ‘ 


oooo 


3 ,*  INNER  = *»rl0.7t/‘  -ti«:u, 

4 /»  • ,T34, ‘NUMBER  = ‘,15,//) 
il5  FOPMATI ‘0*»T20,*Q( 12 


A/ 1 r^uUr**' c r\ciN  I Ji  HC  nuuc 

iREAL  = * jFi0.5t/‘0* t 


C 

C 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


l T20, 'Q(  S I2t*  f • »I2,‘ HMAG  = ',F10.5t/) 

136  FORMAT! 'I') 

END 

THIS  IS  SUBROUTINE  INPUT 
SUBROUTINE  INPUT (lER) 

SUBROUTINE  INPUT  READS  INPUT  VARIABLES,  AND  DEFINES 
SEVERAL  CONSTANTS  USED  THROUGHOUT  THE  PROGRAM. 

COMMON/ BLCKl/  FSTRMN, DELTAS ,DX, DH,DSTSTR , RO ,R I ,REDFRQ 
COMMON/ 3LCK2/  A1 , A2 , A3 , A4 , R N,RM , RM 1 , F ACT  1 , FACT2 
CCMM0N/BLCK3/  NGRDFN , N PTS, LCCUNT , IHAVEP , N, M , Ml , I PR INT 
NAMEL I ST/NAM  1/  FSTRMN, REDFRQ ,RO, RI ,NGRDF N ,N ,M,Ml , 

1 IPRINT 

WRITE(6,80) 

READ(5,NAM1) 

WRIT£(6,NAM1 ) 

FMANGL  = A RSINd  ./FSTRMN) 

DH  = (RO  - Rn/FL0AT(N6R0FN] 

IF  ( DH  .LT.  0.  ) GO  TO  1 

DELTAS  = OH=^FSTRMN 
DX  = 0H/TAN( FMANGL) 

DSTSTR  = 2.*OX 
NPTS  = 1. /'DSTSTR  + I 

IF  ( NPTS  .GT.  400  ) 60  TO  2 

BETA  = SQRT( FSTRMN^FSTRMN  ~ 1.) 

RN  = FLCAT(N) 

RP  = FLOAT (M) 

RPl  = FL0AT(M1 ) 

FACTl  = RE0FRQ=5=REDFRQ 
FACT2  = RN'-^RN/IFSTRMN^FSTRMN) 

COMPUTE  CONSTANTS  TO  BE  USED  IN  THE  COMPUTATIONAL 
MOLECULES. 

A1  = 0.25*0ELTAS/FSTRMN 
A2  = 0.5=t<FSTRMN/BETA 
A3  = REDFRC*DELTAS*A2 
A4  = 0.5-FSTRMN 
RETURN 

1 lER  = 1 

WRITE(6,100) 

GO  TO  3 


C 

C 

C 

c 

c 

c 


,:ABNORMAL  TERMINATION:  INPUT  RADII  •, 
,35X,‘ ABNORMAL  TERMINATION:  MUST  INCREASE 


2 lER  = 1 
WPITE(6,110) 

3 RETURN 
80  FORMAT! ‘1* ) 

100  FORMAT!//*  • ,35X,‘ 

1 'INCORRECT*,//) 

110  FORMAT!//*  * ' 

1 *E  THE  SIZE  OF  VECTORS  XPT,  CPR,  AND  CPI*,//) 
END 

THIS  IS  SUBROUTINE  INITAL 
SUBROUTINE  INITAL 

SUBROUTINE  INITAL  INITIALIZES  ALL  FLOW  FIELD 
QUANTITIES. 

COMPLEX  DPHIOR,DPHIDX 


M 


no 


oooo  o onoo  no  o on  noon 


♦ CkH^OO) 


CCMPLEX  PKK^OO)  i DPHISH4001  ,DPHIS2(400) 
dimension  X(400) »R(400) ,XPT(400»  tCPRl400  i 
COMMON/ BLCKl/  FSTRMNtOrLT &S t DX, DH iDSTST R . RO ,RI , REDFRQ 
Al,A2,A3tA4,RN,RM,RMl,FACTl,FACT2 
NGRDFN«NPTS,LCOUNT, IHAVEPtNtMtMl, IPRINT 
PHI ,0PHIS1,DPHIS2 
X,R»XPT,CPR,CPI 


COMMON/ BLCK2/ 
COMMON/ BLCK3/ 
CCMM0N/BLCK4/ 
COMMON/ BLCK5/ 


DATA  PI/3.141593/ 

xm  = 0. 

R(l)  = RQ 
XPT(l)  = 0. 

ALONG  THE  INITIAL  RIGHT  RUNNING  CHARACTERISTIC  PHI  AND 
0PHIS2  ARE  ZERO. 

PHim  = {0.,0.) 

0PHIS2(1)  = (0.,0.) 


AA=.2/FSTRMN 

OPHISU  U = CMPLX(AA,0.) 

OPHIOR  = A4*0PHISi(l ) 

OPHIOX  = A2^0PHIS1(1 J 
CPR(l)  = -2.=!‘REAL(0PHIDX} 

CPKU  = -2.*AIMAG(0PHIDXI 

IF  ( IPRINT  .NE.  1 ) RETURN 
WRITE(6,90) 

WRlTE(6t95) 

WRITE (6, 100)  Xm,OPHISl(l)  »DPHIS2m  ,PHI(  1), 

I R(l) , OPHIOR, DPHIOX 
WRITEI6,110)  CPR{ 1),CPI m 
RETURN 

90  FORMAT  (‘1‘  ,T  10 , « X S T22 , ♦ 0<  PH  I)/D(  SU  » , T46, 

1 ’0(PHI)/0(S2)* ,T70,'PHI‘ ) 

95  FORMAT!*-  • ,T 10 , • R • ,T23, * D( PHI ) /DR* ,T47, » D( PHI ) /OX* , / ) 
100  FORMATOH  I , 2X,  F8  .454X,  3 ( F 1 2 . 5 , F9 . 5 ,3X  ) ,/ * »,4X,F8.4, 
1 4X,3(F12.5»F9.5,3Xn 

llO  FORMAT! ‘0* ,10X,* INITIAL  VALUE:  INNER  SURFACE  OF*, 

1 * OUTER  CYLINDER*, /•  * ,30X, ’CPIREAL)  = *,F8.4,10X, 

2 *CP{  IMAG)  = * ,F8.4, //) 

END 

THIS  IS  SUBROUTINE  COMPXR 
SUBROUTINE  COMPXR 

SUBROUTINE  COMPXR  COMPUTES  THE  (X,R)  COORDINATES  OF  A 
GRID  POINT. 

01  MENS  I CN  X(400 ) ,R (400 ) ,XPT(400 ) ,CPR( 400 ) ,CPI(400) 
CCMMON/BLCKl/  FSTRMN, DELTAS, OX, OH, DSTSTR, RO,RI , REDFRQ 
CGMM0N/BLCK3/  NGRDFN,NPTS, LCOUNT , IHAVEP , N,M , HI , I PRINT 
CGMMQN/BLCK5/  X, R, XPT ,CPR, CP  I 


NEW  RIGHT  RUNNING  MACH 


IF  THIS  IS  THE  FIRST  POINT 
LINE,  GO  TO  1. 
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i 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 


SUBROUTINE  MACHLN 

SUBROUTINE  MACHLN  COMPUTES  THE  VALUES  OF  THE  FLOW 
FIELD  QUANTITIES  ALONG  THE  INITIAL  RIGHT-RUNNING  MACH 
LINE  WHERE  PHI  AND  DPHIS2  ARE  ZERO. 

COMPLEX  A,C 

COMPLEX  DPHIDRtOPHIOX 

COMPLEX  PHK400)  tOPHISl(400)  ,DPHIS2(400» 

DIMENSION  X(400) , R (400 1 »XPT( 400  ) ,CPR( AOO ) t CP  I ( AOO » 
CCMMGN/8LCK1/  FSTRMN,OELTAS ♦ DXt DHtDSTSTR , RO , RI t REDF RQ 
All A2t A3,A4,RN,RM,RM1 ,FACTl, FACT2 
NGRDFN,NPTStLCOUNT,IHAVEP»N» H,M1, IPRINT 
PHI ♦DPHIS1,0PH1S2 
X,R^XPT,CPR,CPI 


COMMON/ BLCK2/ 
COMMON/ 3LCK3/ 
COMMON/ BLCK4/ 
COMMON/ BLCK5/ 


WITH  THE  RIGHT- 


C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 


I = IHAVEP  *■  1 
J = IHAVEP 

CALCULATE  THE  COEFFICIENTS  ASSOCIATED 
RUNNING  MACH  LINE. 

LINE. 

AAl  = 1.  - A1/R(I ) 

A = CMPLXIAAl, A3) 

Cl  = 1.  A1/R(J) 

C = DPHIS1{J)^=CMPLX(C1,-A3) 

PHKI)  = (0.,0.) 

0PHIS2(I)  = (0.,0.) 

DPHISK  I ) = C/A 

DPHIDR  = A4=!‘0PHISl{n 
OPHIOX  = A2«0PHlSl(n 

IF  ( IPRINT  ,NE.  I ) RETURN 

WRITE  (6f  100)  X(I)  ,0PHISl(I),DPHIS2(n,PHI(I  ). 

1 R( I ) tDPHIDRtDPHIDX 
WRITE(6ail) 

RETURN 

100  FORMATOH  M , 2X,  F 8 .4 , 4X»3  ( F 12. 5 » F9 . 5 1 3X ) » / ' »,4X,F8.4, 
1 4X»3(F12.5,F9.5»3X) ) 

111  FORMAT! ’O'  ) 

END 

THIS  IS  SUBROUTINE  RADI 

THIS  IS  SUBROUTINE  RADI 
SUBROUTINE  RADI 

SUBROUTINE  RADI  COMPUTES  THE  FLOW  FIELD  QUANTITIES  AT 
THE  SURFACE  OF  THE  INNER  RADIUS  WHERE  THE  3CUNDARY 
CONDITION  IS  DEPENDENT  ON  N , THE  CIRCUMFERENTIAL 
MODE  NUMBER.  IF  N IS: 

EVEN  — DPHIDR  = 0,  AND  THUS  DPHISl  = DPHIS2. 

ODD  — CP  = 0. 

COMPLEX  A,B,C»TANRf DPHIORtDPHI DX,DELTA»GAMMAtIK 
COMPLEX  PHK400)  , OPHI  SI  (400  ) ,OPHI  S2(  400  ) 

DIMENSION  X(400i t R (400) ♦ XPT( 400 ) ,CPR(400 ) , CPI ( 400 ) 
CCMM0N/8LCKI/  FSTRMN, DELTAS, DX, DH,DSTSTR, RO^RI , REDFRC 
Al ,A2,A3,A4,RN,RM,RM1,FACT1,FACT2 
NGROFN,NPTS,LCOUNT, IHAVEP, N,M, Ml, IPRINT 
PHI ,DPHISl,0PHIS2 
X,R,XPT,CPR,CPI 


CGMM0N/BLCK2/ 

CGMM0N/BLCK3/ 

CCMM0N/8LCK4/ 

C0MM0N/BLCK5/ 


C 

C 

C 

C 


I = IHAVEP  f 1 
J = IHAVEP 

CALCULATE  THE  COEFFICIENTS  ASSOCIATED  WITH  THE  RIGHT- 
RUNNING MACH  LINE. 
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‘ V _,*-«’'^*»  ff '♦“■'M  W*-  ■*■?'-• 


«»*»••  #••-**  »/-».»»*»^  ««Ar%i^A^>r>  V*/  »*A<»  r-^iiiir^  •ir/^  rwr* 

inc  ruLLUWiiNVa  VHr^iHoucf  ii\»  MHO  ruuiwL»  lu  dc 
UNDEFINED  IN  THE  ORIGINAL  REFERENCE  12  PROGRAM,  AND 
WAS  THE  REASON  FOR  A NUMERICAL  INSTABILITY  AFTER  THE 
FIRST  REFLECTION. 

IK=CMPLX(0.0,RE0FRQ) 

AAl  = 1.  - Al/RU) 

A = CHPLXI  AAl,  A3I 

B1  = Al/Rin  - 0.25=«‘OELTAS=i'DELTAS*(FACTl  - 
I FACT2/(R(  D’S'RI  IJ)) 

B = CMP  LX (Bl, A3) 

Cl  = 1.  + Al/RU) 

C2  = -Al/RU:  *■  0.25*0£LTAS*DELTAS*(FACTl  - 
1 FACT2/I  R(  I)«R(I ))  ) 

C3  = (2.«FACT1  - FACT2/(R( J)«R( J))  - 
1 FACT2/(R(  I)=!=Rm))*OELTAS*0.5 
C = OPHISl  U)«CMPLX(C1,-A3)  + DPHIS2  U ) =<'CMPLX(C2,-A3  ) 
1 + PHI(J)*CMPLXIC3,0.) 


H00{N,2)  .EQ. 


GO  TO 


APPLY  THE  BOUNDARY  CONDITION  FOR  ODD  VALUES  OF  N. 

CP  = (IK) PHI  + DPHIDX  = 0.  USE  THIS  B.C.,  THE  FINITE 
DIFFERENCE  EQN.  FOR  OPHIDX,  AND  THE  EtiN.  FOR  THE  RIGHT 
RUNNING  MACHLINE  TO  SOLVE  FOR  DPHISl  AND  DPHIS2. 

DELTA  = A2  + 0.5*IK*DELTAS 

GAMMA  = IK*(  PHKJ)  <•  0 . S’^'DELTAS^DPH  I S2  ( J)  ) 

% 

SOLVING  FOR  DPHISl  AND  DPHIS2 

OPHISUI)  = { C + 8*GAHMA/DELTA  )/(  A - B’l'AZ/DELTA  ) 
DPHIS2(I)  = ( C - A’^'DPHISim  )/B 

DPHIDX  = A2>!'-(  DPHISKI)  * DPHISZd)  ) 

DPHIOR  = A4=!^(  DPHISKI)  - DPHIS2(I)  ) 

GC  TO  2 

APPLY  THE  BOUNDARY  CONDITION  FOR  EVEN  VALUES  OF  N. 


I DPHISKI)  = C/(A  + B) 
0PHIS2(  I)  = DPHISKI  ) 
DPHIDR  = (0,,0.) 

DPHIDX  = A2^(DPHISKI) 


DPHIS21 I )) 


2 PHKI)  = PHI(J)  + 0.i>*(DPHIS2(  J)  + DPHI  S2(  I } ) ^DELTAS 
IF  ( IPRINT  .NE.  1 ) RETURN 

WRITE  (6, 100)  X(  I ),  DPHISKI  ) ,DPHIS2(  I ) , PHI  I I ) , 

1 R(I ) tDPHIDR, DPHIDX 
RETURN 

100  FORMATOH  RI , 2X,  F8.4,4X,3(  F 1 2. 5 , F9 .5 ,3X  ) , / • •,4X,F8.4, 
1 4X,3(F12.5,F9.5,3X) ) 

END 

THIS  IS  SUBROUTINE  RAD2 
SUBROUTINE  RA02 

SUBROUTINE  RA02  COMPUTES  THE  FLOW  FIELD  QUANTITIES  AT 
THE  INNER  SURFACE  OF  THE  OUTER  CYLINDER,  WHERE  THE 
FLOW  TANGENCY  CONDITION  PRESCRIBES  DPHIDR. 

COMPLEX  D,E,F,TANR 
COMPLEX  DPHIOR, DPHIDX 

COMPLEX  PHI  (400)  ,DPHI  SK400)  ,0PHIS2(400  ) 

DIMENSION  X(400) ,R (400) , XPT( 400 ) ,CPR( 400 ) ,CPI ( 400 ) 
COMMON/ BLCKl/  FS7RMN, DELTAS, DX, DH,0STSTR ,R0 ,RI , REDFRQ 
COMMON/ BLCK2/  A1 , A2, A3 , A4, RN, RM, RMl , FACT  1 , F ACT2 
CCMMON/BLCK3/  N6RDFN,NPTS , LCOUNT , IHAVEP ,N, M ,M1 , 1 PRINT 


OOOOOO 


oooo  oo  n ooo  ooo  o ooooooooo  ooooo 


CCMM0N/BLCK4/  PHI tDPHI SI , DPHIS2 
C0MM0N/BLCK5/  Xt R » XPT ,CPR t C P I 
DATA  PI/3.141593/ 

I = IMAVEP  4-  1 
J = IHAVEP  2 
K = LCOUNT 
ABG=RM*PI-^X(  U 


CALCULATE  THE  CCEFFICIENTS  ASSOCIATED  WITH  THE  LEFT 
RUNNING  MACH  LINE. 

D1  = Al/Rin  + 0.25*DELTAS=i=DELTAS’S'(FACTl  - FACT2/(R(I) 
1 =«=R(I))J 

D = CMPLX(01t-A3J 
El  = -1.  - A1/R(I) 

E = CMPLXIEl ,-A3) 

FI  = -A1/R(J}  - 0.25*DELTAS=«'DELTAS*IFACT1  - FACT2/ 

1 (R(I)^R(nn 
F2  = -i.  + A1/R(J) 

F3  = (2.=!'FACri  - FACT2/(R(  I)*R{  ! J ) - FACT2/(R(J)* 

I R( J) J)^0ELTAS^0.5 

F = DPHISK  J)*CMPLX{F1,A3)  + DPHI  S2(  J1  «CMPLX(  F2,  A31  • 

1 PHI  (J)’5'CMPLX(F3,0.) 


CALCULATE  THE  FLOW  TANGENCY  CONDITION  - 
BOUNDARY  CONDITIONS  FOR  PANEL  FLUTTER 


TANR 


THE  FOLLOWING  THREE  CARDS  WERE  MODIFIED  FROM  THE 
ORIGINAL  REFERENCE  12  PROGRAM. 

ALPH'\=.  1 
TANR1=ALPHA 

TANR2  = REDFR0*Xm*  ALPHA 
TANR  = CM0LX(TANR1»TANR2) 


DFHISim  = ( F <-  TANR*E/A4  )/(  D + 
0PHIS2(I)  = ! F « TANR=<^D/A4  )/(  D + 
DPHIDR  = A4^(0PHISl( I ) - DPHIS2(I)) 
DPHIDX  = A2*(0PHIS1( I)  + DPHIS2(I>) 


INTEGRATE  ALONG  THE  LEFT  RUNNING  MACH  LINE  TO  FIND  PHI 
PHKIl  = PHKJ)  + 0.5*(DPHIS1(  J)  + DPHI  SK  I ) )=«'DELTAS 
CALCULATE  THE  PRESSURE  COEFFICIENT  AT  THIS  GRID  POINT. 

XPT(K)  = xm 

CPP.(K)  = -2.*(-REDFRQ*AIMAG(PHI(  I n + REAL!  CPHIDX)  ) 
CPKK)  = -2.*(HE0F'  Q=!'REAL(PHn  I)  ) <■  AIMAGI  OPHIOXU 

IF  ( IPRINT  .NE.  1 ) RETURN 

WRITE  (6»  100)  XU),DPHISim  ,DPHIS2mtPHim, 

I R( n ,DPHIDR»DPHIDXtTANR 
WRITH(6»110)  CPR(K)tCPI IK) 

RETURN 

100  FORMATOH  R0»  2X  t F8.4,  4X,  3(  FI  2 . 5 » F9.5 ,3X  ) ,/ • *,4X,F8.4, 
1 4X,3(F12.5,F9.5,3X) ) 

no  FCRHAH  »0*  ,10X,*  INNER  SURFACE  OF  OUTER  CYLINDER',/* 

1 30X, 'CP(REAL)  = ' ,F8.4tlOX, 'CP( IMAGI  = »,F8.4,//) 

END 

THIS  IS  SUBROUTINE  6ENFPT 
SLbROUTINE  GENFPT(IER) 

SLBROUTINE  GENFPT  COMPUTES  THE  FLOW  FIELD  QUANTITIES 
AT  A GENERAL  FLOW  FIELD  POINT. 


CCPPLEX  PHIL 


oooo 


oooo  ooooo  o nooo  oo 


CCNPLEX  A»B,C,OtE,F,OENOM 
COMPLEX  DPHIDP jDPHIDX 

COMPLEX  PHI (400) , DPHI SI ( 400 ) ,DPH I S2( 400) 

DIMENSION  X(400»  , R (400  ) » XPT  ( 400  ) ,CPr<(400  ) , C PH  400  i 
COMMON/ BLCKl/  FSTRHN, DELTAS , DX, OHtOSTST R, RO, RI , REDFRC 
Alt A2 » A3,A4, R Nf RM, RMIt FACT 1,FACT2 
NGRDFNtNPTStLCOUNTt IHAVE?,NtM,Ml»IPRINT 
PHI,DPHISl,DPHIS2 
X,RtXPT,CPR,CPI 


C0MMQN/BLCK2/ 
CCMM0N/BLCK3/ 
COMMON/ BLCK4/ 
COMMON/ BLCK5/ 


I = IHAVEP  + 1 
J = IHAVEP 
K = IHAVEP  * 2 

COMPUTE  THE  COEFFICIENTS  OF 
FACT3  = FACT2/(R(  I)*R(in 
FACT4  = FACT2/(R(K)*R{KM 
FACTS  = FACT2/(R(  J)=«'R(  J) ) 


DPHISi  AND  DPHIS2. 


AtB  AND  C ARE  ASSOCIATED  WITH  THE  RIGHT  RUNNING  MACH 
LINE,  WHILE  D,E  AND  F ARE  ASSOCIATED  WITH  THE  LEFT. 

AAl  = 1 . - Al/R( I ) 

A = CMPLXIAAl, A3) 

B1  = -AAl  + 1.  - 0.25=«'0ELTAS*0ELTAS«(FACT1  - FACT3) 

E = CMPLX(B1,A3) 

Cl  = 1.  + A1/P(J) 

C2  = -Cl  4-  1.  + 0.25*DELTAS*DELTAS^(FACT1  - FACT3) 

C3  = (2.=f'FACTl  - FACT3  - FACTS  I *DELTAS‘^0  .S 
C = DPHISK  J)*LMPLX{C1,-A3)  + DPHI S2  ( J)  =i'rMPLX(C2 , -A3 ) 
1 + PHK  JJ^i'CMPLXICSfO.) 

D1  = A1/R(I)  f 0 .25^DELTAS*0ELTAS»!‘(FACT1  - FACT3) 

0 = CMPLX(Dl,-A3 J 
El  = -1.  - A1/R(I) 

E = CMPLX(E1,-A3) 

FI  = -A1/R(K)  - 0.25’!'DELTAS*DELTAS#(FACT1  - FACTS) 

F2  = -1.-  <■  A1/R(K) 

F3  = <2.*FACT1  - FACT3  - FACT4)^DELTAS=i=0.5 
F = 0PHIS1(K)^=CMPLX(F1,  A3)  * U PH  I S2  ( K ) -SC  MPLX(  F2,  A 3 ) - 
1 PHI(K)*CMPLX(F3,0.) 


USE  CRAMERS  RULE  TO  SOLVE  FOR 
PRCVIDED  THAT  THE  DETERMINANT 
SINGULAR, 


DPHISI  AND  DPHi; 
OF  COEFFICIENTS 


is 


NGN- 


DENGM  = A*E  - D*B 

IF  ( CABS(  DENOM  ) .LT.  l.E-05  ) GO  TO  1 
DPHISKI)  = (C*E  - B*F)/DENOM 
DPHIS2tI)  = (A=^F  - C*D  )/DENOM 
DFHIDR  = AA^J^IOPHISHI)  - DPHIS2U)) 
OPHIDX  = A2*(DPHISl(  I ) + DPHIS2(D) 


INTEGRATE 

PHI. 


ALONG  THE  RIGHT  RUNNING  MACH  LINE  TO  FIND 


PHI(I)  = PHHJ)  * 0,5*(DPHIS2{  I ) + DPHI  S2(  J ) )=^<DELTAS 
PHIL  = PHI(K)  + 0,5=!'(DPHIS1(K)  + DPHISiU  ) ) ^DELTAS 

IF  ( IPRINT.NE.  1 ) RETURN 

WRITE  (6, 1 00)  X(  I ) ,DPHIS1U  ) ,DPHIS2(  I ) , PH  I ( I ), 

I R(I J ,OPHIDR,DPHIOX, PHIL 
WRITE<6,111) 

RETURN 

1 lER  = 1 

WRITE(6,110)  DENOM 
RETURN 

100  F0RMAT13H  G , 2X, F 8 .4, 4X,3 ( F 12. 5 , F9.5 ,3X ) ,/ • *,4X,F8.4 
4X,3(F12.5,F9.S,3X) ) 


110  FORMAT! «0* ,///,»  », 35X ABNORMAL  T 

OF  COEFFICIENTS  IS 


'DETERMINANT 


"NATION: 


ngular»,/'5* 


sox 
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WJJ  w«4.*i  W..IM1H14UP 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


2 'DENOMINATOR  = »,E12.5) 
111  FORMAT! 'O' ) 

END 

THIS 

SUBROUTINE  QSF 
SUBROUTINE  QSFl H,Y»Z »NDIM) 
CALL  QSF  (H,Y,Z,NDIM) 
DESCRIPTION  “■  " ■ 

H 


C 

C 


C 

C 

C 

C 


IS  SUBROUTINE  QSF 


Y 

Z 


NDIM 


nP  parAm'^tfr'J 

THE  INCREMENT  OF  ARGUMENT  VALUES. 

THE  INPUT  VECTOR  OF  FUNCTION  VALUES. 
THE  RESULTING  VECTOR  OF  INTEGRAL 
VALUES.  Z MAY  BE  IDENTICAL  WITH  Y. 
IDENTICAL  WITH  Y. 

THE  DIMENSION  OF  VECTORS  Y AND  Z. 


REMARKS 

NO  ACTION 


IN  CASE  NDIM  LESS  THAN  3. 


DIMENSION  Y(1),Z(1) 

HT=,3333333*H 
IF(NDIM-5) 7»8, 1 

NDIM  IS  GREATER  THAN  5.  PREPARATIONS  OF  INTEGRATION 
LCCP 

1 SUMi=Y(2)+Y(2) 

SUMl=SUMlf SUMl 

SUM1=K.  *(YmfSUHl+Y(3n 
AUX1^Y(AHY(4) 

AUXl=AUXl+AUXi 

AUX1=SUM1+HT*(Y(3)4-AUXH-Y<5  ) ) 

AUX2=HT=i'{  V(1  )l-3.  a75#(Y(2)<-Y(5)  )+2.625*{Y(3)  +YI4)  ) + 

I Y(6)) 

SUM2sY(5M-Y{5) 

5UM2=SUM2+SUM2 

SUM2=AUX2-HTt { Y ( 4} +SUM2+Y( 6 ) ) 

Z(1)=0. 

AUX=Y(3)^-Y<3) 

AUX=AUX+AUX 

Z{2)=SUM2-HT=5<(Y(  2H-AUX+-Y{4)  ) 

Z(31=SUM1 

ZI4)=SUM2 

IF(NDIM-6)5^5,2 

INTEGRATION  LOOP 

2 DC  4 I = 7tNDIM,2 
SUM1=AUXI 
SUM2=AUX? 

AUX1=Y(  I-n^-Y{  I-l) 

AUXi=^AUXlfAUXl 

AUXl=SUMi+HT=i‘(Y(  1-2  ) «-AUXli-Y(  I ) ) 

Z(  I-2)=SUM1 
IF( I-NDIM)3,6.6 

3 AUX2=Ya)  + YUi 
AUX2=AUX2+AUX2 

AUX2=SUM2+HT=l*(Y(  I-l  ) + AUX2+Y(  I + l ) ) 

4 Z(I-n  = SUM2 

5 Z(NOIM-l J=AUX1 
Z(i\DIM)=AUX2 
RETURN 

6 Z{NOIM-l)=SUM2 
Z(i\'DIM)^A'JXl 
RETURN 

END  OF  INTEGRATION  LOOP 

7 IF(NDIM-3) 12,11, 8 

NDIM  1$  EQUAL  TO  4 OR  5 

8 SUM2=1.12  5=!'HT=i'(Y{l  UY(2)+Y{2)<-Y(  2 > +Y( 3 ) +YJ  3 H-Y(  3)  *■ 


1 Y(4)) 

SUMl=Y{2M-Y(2» 

SUf^l  = SUMl  + SUMl 
SUMl=HT=i'{Y(  n+SUMl+Y<  3)  ) 
Z(1)=0. 

AUXl=Y{3)<-Y(3) 

AUXl=AUXH-AUXl 

Z(2)  = SUM2-HT=J=(Y(  2M-AUX1+Y(4) ) 
IF(NDIM-5) 10,9,9 
9 AUX1=Y(4)*-Y(4) 

AUXl=AUXl+AUXl 

Z(5)  = SUM1  + HT-(Y(3)+AUXU-Y(5)) 
10  Z(3)=SUM1 
Z(4)^SUM2 
RETURN 


NDIM  IS  EQUAL  TU  3 

U S‘JK1=HT*(  1.25*Y(  1H-Y(2)+Y(2)-.25*Y{3)  ) 
SUM2=Y( 2)+Y( 2) 

SUM2=SUM2+-SUH2 

Z(3)  = HT=«'(Ya  J+SUM2+Y(3H 

Z{1)=0. 

Z(2)=SUM1 
12  RETURN 
END 


1^ 

'm 

\%i 


The  outer  flowfield  program,  from  reference  13,  consists 
of  5 subroutines  and  a main  program  to  manipulate  the 
subroutines.  The  action  of  each  particular  subroutine  is 
indicated  in  the  program  listing.  Inputs  required  to  the 
program  must  be  on  tvjo  cards,  and  are  as  follows: 

DATE 

HACK  Hach  number. 

RADIUS  Radius  to  length  ratio. 

K Reduced  frequency. 

HFRSQ  Axial  mode  number. 

NR  The  r in  Qmr,  the  generalised  aerodynamic  force. 

NFREQ  Circumferential  mode  number. 

FINGRD  Grid  fineness,  essentially  the  number  of  grid 
points  along  the  cylinder  surface. 

The  run  date  is  typed  onto  the  first  card.  The  second  card 
contains  the  above  parameters  in  the  following  format: 
F20.8,  2F10.3,  415.  3f  the  above  parameters,  in  the 

modified  program,  KFRFQ,  NR,  and  NFREQ  must  be  set  to  1,  as 
the  case  under  consideration  is  not  panel  flutter  but  a 
cylinder  at  angle  of  attack.  Other  parameters  are 
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